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ABSTRACT 



Over the last few years, numerical models of the behav- 
ior of solar magnetic flux tubes have gone from using 
methods that were essentially one-dimensional (i.e. the 
thin flux tube approximation), over more or less idealized 
two-dimensional simulations, to becoming ever more re- 
alistic three-dimensional case studies. Along the way a 
lot of new knowledge has been picked up as to the e.g. the 
likely topology of the flux tubes, and the instabilities that 
they are subjected to etc. Within the context of what one 
could call the "flux tube solar dynamo paradigm," I will 
discuss recent results of efforts to study buoyant magnetic 
flux tubes ascending from deep below the photosphere, 
before they emerge in active regions and interact with the 
field in the overlying atmosphere (cf. the contributions by 
Boris Gudiksen and Ake Nordlund): i.e. I am not address- 
ing the flux tubes associated with magnetic bright points, 
which possibly are generated by a small-scale dynamo 
operating in the solar photosphere (cf. the contribution 
by Bob Stein). The presented efforts are numerical MHD 
simulations of twisted flux ropes and loops, interacting 
with rotation and convection. Ultimately the magnetic 
surface signatures of these simulations, when compared 
to observations, constraints the dynamo processes that are 
responsible for the generation of the flux ropes in the first 
place. Along with these new results several questions pop 
up (both old and new ones), regarding the nature of flux 
tubes and consequently of the solar dynamo. 

Key words: Sun; magnetic fields; flux tubes; interior; 
granulation. 



1. INTRODUCTION 



Buoyant magnetic flux tubes are an essential part of the 
framework of the current theories of dynamo action in 
both the Sun and solar-like stars: it is believed that when 
formed near the bottom of the convection zone (CZ), 
by a combination of rotation and turbulent convection, 
toroidal flux tubes buoyantly ascend in the form of tubu- 
lar il-shaped loops. Rising under the influence of rota- 



tional forces, they finally emerge after a few months as 
slightly asymmetric and tilted bipolar magnetic regions at 
the surface. Many models of buoyant magnetic flu x tubes 
are b ased on the thin flux tube approximation (Spruit, 
1981[ ) that treats the tubes as strings, much thinner than 



e.g. the local pressure scale height, moving subjected the 
Coriolis and drag forces. This essentially 1-d approxi- 
mation is consistent with the observations when used to 
study the latitudes of emergence, tilt angles, and the tilt- 
scatter of bipolar magnetic regions on the Sun, but only if 
the initial field strength of the tubes are of the order of 10 
times the c onvection equipartition value near the bottorn 
of the CZ (t)'Silva & Choudhuril , |1993| ; [Fan et al.j [l994 



Caligari et al, 1995). A major problem in dynamo the- 



ory is to understand how the field strength can become so 
high, corresponding to an energy density 100 times larger 
than that of the available kinetic energy. 

However, as buoyant flux tubes rise they expand and the 
assumption that they are thin breaks down some 20 Mm 
below the solar surface. Traditionally, the latter fact is 
taken as the main reason for "going into higher dimen- 
sions", i.e. for submitting to 2-d (actually 2.5-d) and 3- 
d models. So far, a lot of questions still remain unan- 
swered, e.g. whether the quasi-steady state topology that 
the flux ropes reach in the later phase of their rise in 2- 
d simulations (e.g. Emonet & Moreno-Insertis 1998 and 
Dorch et al. 1999) is stable towards perturbations from 
the surroundings, and whether the results found for 3-d 
flux ropes moving in a 1-d average static stratification, 
at all are valid in the more realistic case. Hence when 
attempting to review the status of buoyant magnetic flux 
tube models, there are several questions (Q's) that one 
may ask. Below are a few examples: 
Ql: In general the tubes may be twisted, but how much 
twist is needed and warranted? 

Q2: How does the tube's twist evolve as they ascend: do 

they e.g. kink due to an increasing degree of twist? 

Q3: Are there other instabilities besides the magnetic 

Rayleigh-Taylor (R-T) and kink instabilities? 

Q4: How do the tubes interact with the flows within the 

CZ and at the surface? 

Q5: What happens to the less buoyant magnetic subsur- 
face structure as the tubes rise (the wake)? 
Q6: What happens when the ropes become thick, i.e. 



large comparable to the local pressure scale height? 
Q7: What happens at emergence? Does the field topol- 
ogy change (how)? 

Q8: How does the twist arise? What is the appropriate 
initial condition? 

In the following I will try to answer some of the above 
Q's by reviewing some of the main results of 2-d and 3-d 
models of buoyant magnetic flux tubes with emphasis on 
the most recent results from the turn of this millennium. 



2. REFERENCE MODELS 



For my discussion of the various aspects of flux tube 
models in the following, I have made a number of 2-d 
and 3-d simulations for reference. These are simulations 
for several values of the characteristic field line pitch an- 
gle V'R- 

The 2-d reference simulations lack convection and the 
flux tubes move in a polytropic atmosphere (dynamic), 
while the 3-d reference models retain full hydrodynamic 
convection with solar-like super-granulation and down- 
flows (see the Appendix at the end of this paper). 




Figure 1. Flux tube disrupted by the magnetic Rayleigh- 
Taylor instabili ty; from a 2-d simulation by E monet & 
Moreno-Insertis ( 199^ : The typical pitch angle was 
set to 7° (left), 2°. 5 (middle) and 0° (right). Shown 
are equipartition curves (top row), field strength (bottom 
row), and flow pattern (tiny arrows). 



3. TWISTED TUBES 



To identify the models I use a terminology where the 
models are referred to as e.g. 3D25 for a 3-d model where 
the flux tube initially is twisted corresponding to a pitch 
angle of ^Pr — 25°. A 2-d model with the same pitch an- 
gle is referred to as 2D25. All the reference models have 
plasma-/3's at the onset of the model equal to (3o = 100. 
The set-up of the reference models are discussed in more 
detail in the Appendix, but the initial twist of the flux 
tubes is given by 



= Boe-(''/^)' and = ar/R 



(1) 



where is the parallel and B^ the transversal compo- 
nent of the magnetic field with respect to the tube's main 
axis. Bo is the amplitude of the field, R the radius, and 
a is the field line pitch parameter. The tubes are twisted, 
with a maximum pitch angle of t/jr = arctan(Q!), hence 
"rope" is a better word to describe the state of the field 
lines than "tube". In the following I set the initial radius 
to Ro = 0.177 Hpo, i.e. the tubes are non-thin. 

The wavelength A of the flux rope is equal to the hor- 
izontal size of the domain A = 3.2 Hpo at the initial 
position of the rope where the pressure scale height is 
Hpo = 78.3 Mm. Then the flux rope is not undular 
Parker-unstable even though the stratification admits this 
instability for ropes longer than a crit i cal w avelength of 
- 12 Hpo ( [Spruit & van Banegooijer] , |1982| ). 



Primarily the behavior of the flux tubes during the initial 
and rise phases are discussed. When the ropes approach 
the upper boundary of the CZ, they enter an "emergence 
phase" slightly beyond the scope of this review (but see 
e.g. the contributions by Gudiksen and Nordlund, ibid). 



2-d simulations of flux tube cross-sections, not relying on 
the assumption that the flux tubes are thin, have shown 
that purely cylindrical tubes are quickly disrupted by 
the magnetic R-T instability, rendering them unlikely to 
reach the surface because thev loose their bu ovancy and 
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cohesion ( Schiisslerl |1979|; Fsinganos 



Mor eno-Insertis, |1998|; porch & Nordlund), 
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et al., |1998| ; [Wissink et alj , pOOO[ ). This apparent disrup- 
tion of the tubes is a result of the too simple topology 
assumed for their magnetic field: rather, the magnetic 
field line tension present in a twisted flux rope suppresses 
the R-T instability and hence prevents the flux structures 
from disintegrating as it has been demonstrated by e.g. 
Emonet & Moreno-Insertis (1998): an expression similar 
to Eq. (1) has been applied in both 2-d and 3-d simula- 
tions and it has been shown that for these cases the R-T 
instability is inhib ited if the degree of twist is suffici ently 
high as shown by ^monet & Moreno-Insertis (1998 ): the 
corresponding critical value of the field line pitch angle 
\E'c is approximately determined by equating the energy 
density of the transversal (or azimuthal) field component 
equal to the ram pressure of the flow relative to the flux 
rope (incidentally, several approaches yields the same ex- 
pression for ^c)- This kind of twist-topology is simple 
but similar to the relaxed state of a flux rope with a mor e 
complicated topology used i ntPorch&Nordlunc^ ( |1998| ). 
Emonet & Moreno-Insertis ( |l99^ finds a critical pitch 
angle on the order of arcsm (-^/R/Hpo), which typically 
is 10° for thin flux ropes. For thicker flux ropes, like the 
ones I discuss here, the value of 4'c increases to approxi- 
mately 25°. 



The boundary of a flux rope can be defined by the 
equipartition curve (the locations where the transversal 



Figure 2. Images showing at time 67.7 TAfor2D30 (left), and at time 35.6 rp^for 3D25 and 3D45 (center and right). 
Shown by contours are the equipartition curves (full curves). Also shown are HWHM-contours ofJi^ (dashed). Only a 
subsection of the full computational box is shown. 



magnetic field is in equipartition with the kinetic energy 
relative to the rise of the rop e). Fig shows equiparti 



t ion cu rves for 2-d models by Emonet & Moreno-Insertis 



( 1998 ), corresponding to three values of the pitch an- 
gle: only the most twisted ropes withstand the exterior 
pressure fluctuations, and the transversal field is strong 
enough to protect the buoyant core of the flux rope. 

The results of my 2-d reference simulations agree per- 
fectly with the 2-d studies mentioned above: the ropes 
are disrupted by the R-T instability unless their pitch an- 
gle exceeds a critical value. As the ropes rise and expand 
they enter a "terminal rise phase" where they ascent with 
a constant speed — while oscillating due to their differen- 
tial buoyancy — ^before reaching the upper computational 
boundary (which is open). In this rise-phase the speed 
is the so-called terminal velocity vt determined by the 
balance between buoyancy and drag. In reference simu- 
lation 2D30— see Fig. | (left)— I get vt « 0.1 dao (wifli 
WAo being the Alfven speed wa at the initial position). 
The largest Mach number reached during the entire rise is 
8.5 X lO^'^, so that the motions are completely sub-sonic. 
The characteristic time scale of the ascent then becomes 
Triso = Hpo/wt ~ 50 TA with TA = RqA'ao- From the 2- 
d reference simulations I find that a minimum pitch angle 
in the range 25° - 30° is required for the ropes to be un- 
af fected by the R-T instability , consi stent with the result 
of Emonet & Moreno-Insertis ( |1998| ). 



Consider a 3-d reference simulation with '^-r = 25° 
(3D25): in that case the pitch is far to low to excite 
the kink instability (see next Section), but we expect to 
re-find the R-T instability in addition to the action of 
the convective flows. Fig |[ (center and right) shows 
cross-sectional averages of the magnetic field in 3D25 
and 3D45 which may be compared to the correspond- 
ing image from 2D30 (note the difference in vertical off- 
set). In the 2-d case, the horizontal and vertical sizes of 
the ropes do not increase equally fast. This effect is pri- 
marily because of the differential buoyancy that squeezes 
the ropes in the vertical direction, and by flux conserva- 
tion, expands them horizontally. The effect is less ap- 
arent in 3-d, where the ropes are more circular (see Fig 
when defined by the HWHM-boundary (i.e. the half- 
width at half-maximum of the parallel field component). 
The equipartition curve, however, loses its 2-d mirror 
symmetry. For a less twisted rope (e.g. 3D25 in Fig. 0), as 



the rope rises and the R-T instability sets in, the flux rope 
disrupts and eventually fills the bulk of the CZ. At late 
stages the magnetic field forms network patches at the 
surface cell boundaries, while the pumping effect (P orch 
& Nordlund, [2001 ) transports the weakest flux down to- 
wards the bottom of the CZ. 



3.1. Rotation and 0-loops 

It has been suggested that the value of the critical degree 
of twist needed to prevent the R-T instability may be un- 
realistically high in the 2-d case, and a smaller twist may 
be sufficient in the case of sinusoidal 3-d m agnet ic flux 
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loops as in the simulations: Abbett et al. 
their anelastic MHD code to study the buoyant rise of 
twisted flux loops and found that the greater the curvature 
of the loop at its apex, the sma ller the critical pitch angle. 
Moreover, Abbett et al. ( 2001 ) found that including solar- 
like rotation and a Coriolis force when solving the MHD 
equations lead to less fragmentation, even of initially un- 
twisted buoyant flux tubes. However, the limit on '^c is 
not known. 



3.2. Sigmoids 



Recently, 3-d simulations of buoyant twisted flux loops 
have confirmed that several of the results found in the 2-d 



(Matsumoto et al. 
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rom the 3-d result 



it is apparent that the apex cross-sections of twisted flux 
loops are well described by the 2-d studies. Furthermore, 
the S-shaped structure of a twisted flux loop as it emerges 
through the upper computational boundary is qualitative 
similar to the sigmoidal structures observed in EUV and 
soft X-ray by the Yohkoh and SoHO satell ites, see Fig. |] 
dCanfield etaL] , |1999| ; ^rUng et al.| , |2000[ ). 



An example of suc h 3-d simulations are provided by M a- 
gara & Longcope ( [200 1[ ) who performed ideal MHD 
simulations of a flux rope emerging from a convection 
zone through a photosphere into an overlying corona: ini- 
tially the flux rope is horizontal and in mechanical equi- 




Figure 3. S-shaped structure of emerging magnetic flux rope from a simulation by ]porch etal. ( 199^ ( left), and a sigmoid 
in X-ray from the Yohkoh satellite ( right). 



librium, but when subjected to the undular Parker insta- 
bihty it is forced to rise and break through the surface, be- 
cause of the buoyancy enhancing downflows alon g the it s 
field lines. A similar study was performed by Fan (2001). 
The result is that first a magnetic bipole is formed at the 
surface by the oppositely directed vertical field lines of 
the looping flux rope. Secondly, upon emergence, a neu- 
tral line and an S-shaped (sigmoidal) structure is formed 
by the flux rope's central (inner) field lines see Fig. Q The 
outer field lines expand and form an S with the opposite 
sign of the inner one, as well as an arcadial structure with 
an almost potential geometry. 



3.3. The flux rope's path 

The naivistic view of a buoyant flux rope ascending 
nicely on a straight vertical (i.e. strictly speaking radial) 
path towards the surface can be dismissed qualitatively; 
in the 3-d case with convection, it is easy to imaging that 
the rope will be influenced by convective flows and be 
pushed around, as we shall see below. However, even 
in 2-d, neglecting flows that is not connected with the 
rope's buoyant rise, a more subtle effect courses the ropes 
to have very complex trajectories when the Reynolds 
number becomes sufficiently large: [Emonet et al. ( 2001| ) 
found a zig-zag motion of the magnetic flux rope caused 
by the imbalance of vorticity in the flux tubes, see Fig. ^ 
Vorticity generated when flux ascend supply a lift force; it 
is by this mechanism that weakly twisted flux ropes loose 
their updrift (negative lift force from two vortex rolls). 
At high viscous Reynolds number the shredding of vor- 
tex rolls in the trailing Karman-like wave of the flux in- 
troduces a vorticity-imbalance, resulting in an alternating 
lift force shuffling the rising rope. Perhaps surprisingly, 
it turns out that the complicated path of the flux rope can 
be described by a single non-dimensional number. 



In the case 3D45, the degree of twist is small enough to 
prevent the onset of the kink instability (the linear ideal 
growth rate vanishes for a = 1, e.g. Fan et al. 1999), yet 
it is large enough to prevent also the onset of the R-T in- 
stability: Fig. mllustrates the effect of these instabilities 
on the shape of the flux ropes. Thus, in 3D45 the rope 
retains its cohesion without distorting its shape by any of 
these two instabilities, and we may focus our attention on 
the effects of the convective flows on the rope. 

The initial position of the ropes in the CZ is significant 



for the subsequent detailed history of their rise: with the 
present flows and location of the flux ropes, most parts 
of the ropes are incidentally located inside or close to 
a convective updraft. Thus, the ascent of the ropes are 
influenced by this: Fig ^ compares 3D45 to 2D45 (i.e. 
to a convection-less reference simulation), and a sim- 
ple analytic flux tube. As the 3-d rope rises, convective 
flows perturb its motion, preventing it from entering a 
rise phase with a constant rise speed, as it indeed does 
in 2D45 (see Fig ^.) The rope remains straight and the 
maximum excursion of its axis, at the end of the simu- 
lation, is ^ 0.04 A. With the chosen super-equipartition 
axial field strength, the main action of the large-scale con- 
vective flows is to push the rope both left and right of the 
central plane (Fig — see also the movie on the CD- 
rom), while the effect of the small-scale downdrafts is to 
locally deform its equipartition boundary. 

The initial location within a general updraft region ex- 
plains why the rise speed of the rope is slightly greater 
than that of the reference simulation, which reaches a ter- 
minal speed of ~ 0.1 vao- The 3-d rope also expands 
more quickly than the rope in the 2-d simulation (see 
Fig. but its rate of expansion is closer to what is ex- 
pected from an adiabatically expanding, non-stretching 
tube with constant flux. 



4. KINKING ROPES 



Fan et al. (1999) report on 3-d non-convective simula- 
tions of buoyant kink-unstable flux ropes in a stratified 
atmosphere, and found that a twisted flux rope becomes 
significantly unstable only when the pitch parameter a 
is well above a critical ideal-MHD value ac = 1, i.e. 
if ip-R, > 45°: in the resistive-MHD case ac increases 
slightly. In reference simulation 3D75 a kink therefore 
develops, with a maximum ideal growth rate of Fkink oc 
r^^ (a^ — 1) (from Linton et al. 1998), yielding a charac- 
teristic time scale of Fj^jJ^,^ w 0.29 ta significantly shorter 
than the rise time derived from the 2-d reference simula- 
tion 2D75, see Fig. |[ Diffusion decreases the growth rate 
further, but the instability still sets in rather quickly, well 
before the rope reaches the surface: the rope has plenty 
of time to kink while rising. 



I do not invoke a perturbation of the rope's mass or en- 
tropy to onset the kink, as do e.g. Abbett et al. (2000) and 
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Figure 5. 2-d simulation of a buoya nt flux rope ri s ing th rough a stratified atmosphere; the rope follows a zigzag trajectory. 
Figure provided by T. Emonet, e.g. Emonet et al. (2001). 




Figure 6. Snapshots of magnetic field isosurfaces at time < = 28 Tp^for the simulations (left to right) 3D25 (weak twist, 
R-T unstable), 3D45 (stable rope) and 3D75 (strong twist, kink unstable rope): the magnetic field isosurfaces outline the 
HWHM surfaces. 



others. Rather I let the rope be perturbed by the convec- 
tive flows. The power of the convection at the initial po- 
sition of the rope is predominantly at long wavelengths 
(100 Mm), and these wavelengths of the unstable kink 
modes dominate the rope's evolution. 

As the rope kinks, it develops a significant curvature 
and a knotty shape with portions dominated by vertical 
field components, see Fig. || (right). Upon emergence 
the bipole rotates and the neutral line between the two 
main polarities are strongly sheared, see Fig. |[ Hence, 
when the rope approaches the surface, it contains adja- 
cent mixed polarities in horizontal cross-sections (see the 
CD-rom movies); this is reminiscent of the appearance of 
(5-type sunspots with different polarities within a common 
umbra, and kinked flux tubes are often quoted as possibly 
being responsible for this phenomenon. 



5. LOOSING FLUX 



[ 



Petrovay & Moreno-Insertis (1997 ) suggested that turbu- 
lent erosion of magnetic flux tubes may take place due 



to the "gnawing" on a flux tube's field lines by turbulent 
convection: the flux tube is eroded by a thin current sheet 
forming spontaneously within a diffusion time. 

Fig ^ shows the 3D45-rope's characteristic radius Rhwhm 
along its axis (with the radius defined as the average 
HWHM). As the rope rises and expands, its magnetic 
field strength also decreases to approximately conserve 
the flux (the decrease found in 3D45 is at a rate close 
but not identical to that expected for a flux-conserving 
tube). The deviation can be attributed to the fact that, 
during its ascent, magnetic flux within the 3-d rope is lost 
to its surroundings. This is illustrated in Fig HQ. (left), 
which shows the total normalized magnetic flux within 
the rope's HWHM-core as a function of time for both 
the 2-d and 3-d ropes (2D45 and 3D45). As the 3-d sim- 
ulation progresses, the total flux-loss from the computa- 
tional domain is only 0.3% — the flux content of the rope, 
however, decreases more rapidly. 

Also shown in Fig |lO[ (right) is the magnetic flux external 
to the rope $e both above and below its center $„ and $; 
respectively. Since the sum $e + is nearly conserved, 
as $i decreases, $e — + must increase by an equal 
amount. However, the distribution of the flux-loss is not 
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Figure 4. Snapshot of emerging twisted flux rope ( dark 
field lines, horizontal velocity field as white arrows, ver- 
tical velocity in gray scale, and vertical magnetic field 
component as contours): The darkest magnetic field line 
is the central axis o f the flux rope. Figure provided by T. 
Magara, e.g. Magara & Longcope (2001). 



symmetric around the rope: more flux is lost to the sur- 
roundings below the rope than above it. 

This asymmetry also exists in 2-d, even though the total 
flux-loss is much smaller in that case. The asymmetry is 
a result of two factors. As the rope rises, the total volume 
above it decreases, while the volume below it increases. 
Furthermore, there is an anti-symmetry of the relative rise 
velocity across the rope: when the rope ascends, there is 
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Figure 7. (a) Height of the flux rope as a function of 
time (3D45, stars and 2D45, diamonds). The straight 
line corresponds to the average speed o/O.l uao in the 
rise phase, (b) Average drift of the 3D45 flux rope in the 
meridional plane. 




Figure 8. Simulation vector magnetogram of an emerging 
kinked rope showing the horizontal ( ar rows) and v ertica l 
magnetic fleld (contours), adopted from Fan et al. ( 1999). 
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Figure 9. Expansion of the flux rope ( 3D45, stars and 
2D45, diamonds) in units of its initial radius. Also plotted 
is an analytical expression (solid line, see text). 



a tendency for flux to be advected towards it near its apex, 
and transported away from it in its wake. The more pro- 
nounced asymmetry in the 3-d case can be attributed to 
the pur nping effect that transpor ts the weak field down- 
wards dPorch & Nordlun4 |200l| ). 



We have defined the flux rope as the core magnetic struc- 
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Figure 10. Left: magnetic flux within the rope $ ; (3D45, 
stars and 2D45, diamonds) and an analytic fit (solid 
curve). Right: the normalized flux outside and above the 
center of the 3D45 rope (stars), and below, (trian- 
gles). The same quantities are shown for 2D45 simula- 
tion (solid and dashed curves respectively). 



ture that lies within the HWHM-boundary. This bound- 
ary is not, however, a contour that moves with the fluid in 
the classical sense: the flux within the latter kind of con- 
tour is naturally conserved (in ideal MHD) and is equal 
to the total flux <f>e + The HWHM-boundary is a con- 
venient way of defining the flux rope and a characteristic 
size Rhwhm, that behaves as it is expected to. The evolu- 
tion of the flux within the rope's core <i>i is determined by 
the integral of Av x B along the rope's boundary, where 
Av is the difference between the fluid velocity and the 
motion of the HWHM-boundary: the average "slip" Av 
is only a small fraction of the rise speed. 

Making the rather crude assumptions that the boundary 
only moves radially relative to the fluid and that the cir- 
cumference of the boundary is circular (which it is cer- 
tainly not, but see Fig. 0), the flux-loss becomes 



'TrRhwhm Aw Be 



(2) 



where Be is the field strength at the center of the flux 
rope. Integrating the equation numerically with the quan- 
tities determined from the simulation 3D45 the result is 
a good fit to the actual flux-loss, see Fig (left)— if 
Av is set to 3 10^'' vaq throughout the time span of the 
simulation except for a short interval of ~ 3 around 
t ~ 30 TA, where At; chan ges sign, as the rop e passes 
from one updraft to another ( porch et~ar , 2001 ). 



Finally, I do not find flux-loss via the type of enhanced 
diffusion in the reference s i mulat ions as proposed by 
^etrovay & Moreno-Insertis ( 1997 ), and rather, the flux- 
loss found here is entirely due to the advection of flux 
away from the core of the flux rope by convective mo- 
tions. The turbulent erosion mechanism requires the tur- 
bulence to be resolved down to scales much below what 
I have implemented here, i.e. my results does not con- 
tradict the existence of the former effect. Most of the 
lost flux ends up in the rope's wake and some is mixed 
back into the upper layers. I speculate that both types of 
flux-loss may take place in the Sun, and as a result, the 



amount of toroidal flux stored near the bottom of the solar 
CZ may currently be underestimated. 



6. DISCUSSION & SUMMARY 



In this review I have tried to shed light on some of the 
"big" questions in the theory of buoyant magnetic flux 
tubes. Below I summarize some of the answers (A's) to 
the questions posed in the Introduction. 



Al: The critical pitch angle is about 30° (for thick ropes), 
in the general 3-d case including (large-scale) convective 
flows, b ut without rotational effects. The result of A b- 
bett et al. ( |200C| ) that the inclusion of rotation lowers the 
critical pitch angle relative to the 2-d result, comes from 
non-convective simulations; i.e. to finally answer this Q 
we need a simulation including both rotation and realistic 
convection. 

A2: I found no evidence that the ropes kink as they rise, 
at least if their initial twist is low enough: one may spec- 
ulate that an initially R-T unstable rope may stabilize due 
to the increase of the twist (which increases only slightly 
in the less symmetric 3-d models). 
A3: There may be other instabilities associated with the 
motion through the uppermost strongly stratified super- 
adiabatic layers, but it is not know from the models pre- 
sented here. 

A4: Besides the shuffling of the flux ropes by the con- 
vective flows, the primary interaction is the cause of a 
flux-loss due to advective erosion. 

A5: The weak magnetic field and the flux lost from the 
rising flux tubes are transported downwards by the flux 
pumping-effect. 

A note on A4: The numerical simulations show that the 
interaction of a buoyant twisted flux rope with stratified 
convection leads to a magnetic flux-loss from the core 
of the rope. During the simulation, the flux rope rises 
96 Mm, and loses about 25% of its original flux con- 
tent. This, ceteris paribus, leads to a small increase in 
the amount of toroidal flux that must be stored at the bot- 
tom of the CZ during the course of the solar cycle: Solar 
toroidal flux ropes rise about 200 Mm before emerging 
as bipolar active regions. One may thus expect them to 
lose even more of their initial flux, which would then be 
pumped back towards the bottom of the CZ. Moreover, 
the relative slip does no t rema in constant throughout the 



rope's rise (Porch et al. 2001 ) 



Of course a lot of Q's still remain to be answered (three 
from my list): the most fundamental problems remaining 
are those of the origin of the twist, and the question of 
how it arises, Q8. This is not addressed by any of the 
models discussed here, but in my view one likely pro- 
cess is the generation of twisted field lines in large-scale 
flux bundles located near the bottom of the convection 
zone, connecting across the solar equator: such flux bun- 
dles would experience a rotating motion since their lower 
parts are located in a region rotating slower than their up- 
permost parts. This rotation would transmit a twist to the 
parts of the flux bundle a slightly higher latitudes, thereby 
possibly giving rise to a twisted toroidal flux system. 
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APPENDIX: SET-UP OF REFERENCE MODELS 



The 2-d reference simulations lack convection and the 
flux tubes move in a polytropic atmosphere calculated by 
using the mass density p{x) and the average gas pres- 
sure P{x) in a 3-d hydrodynamic convection model via 
P = Po{p/ poV^ where Pq and po are the average quanti- 
ties at the initial position of the flux tube. In the 3-d refer- 
ence models the full 3-d convection model is employed. 

The initial set-up of the 3-d models are twofold, consist- 
ing of a snapshot of a solar-like CZ, and of an idealized 
twisted magnetic flux rope. Note that in the 2-d reference 
models, the CZ is polytropic through-out. 

The full MHD-equations are solved on a grid of 150 ver- 
tical times 105^ horizontal grid points, using the compu- 
tational method by Galsgaard and others (Galsgaard & 
Nordlund 1997, and Nordlund, Galsgaard & Stein 1994): 
Typical magnetic Reynolds numbers Rcm in non-smooth 
regions of the domain are of the order of a few hundred, 
while several orders of magnitude higher in smooth large- 
scale regions. 

To model solar-like convection without including all the 
layers up to the actual surface, I use a simple expression 
for an isothermal cooling layer at the upper boundary of 
the model, restricting the effect to a thin layer This layer 
is, however, far below the real boundary of the solar CZ. 
Horizontally the boundaries are periodic. 

The hydrodynamic part of the initial condition is a snap- 
shot from a well developed stage of a numerical "toy 
model" of deep solar-like convection with a gas pressure 
contrast of roughly 2.5 orders of magnitude in the CZ 
alone. The physical size of the computational box is 250 
Mm in the horizontal direction and 313 Mm in the ver- 
tical. A cellular granulation pattern is generated on the 
CZ's surface with a typical length scale of about 50 Mm, 
about twice the canonical size of solar supergranules. The 
typical velocity is about 200 m/s in the narrow downdrafts 
at the surface and slightly less in the larger upwelling re- 
gions, close to what is found for solar supergranulation. 

Initially the entropy in the interior of the tube is set equal 
to that in the external medium. This corresponds to buoy- 
ancy a factor of 1 /j (with 7 = 5 /3) lower than in the 
case of temperature balance where the buoyancy has the 
classical value of 1//3. The initial twist of the flux tubes 
is given by Eq. (jl|). The coordinate system is chosen so 
that X is the vertical coordinate and z the coordinate along 



the axis which initially is parallel to the rope. The posi- 
tion of the rope is described by the set {xc, He), initially 
equal to (a;o,yo)> where the points {xc{z),yc{z)) along 
the z-axis are the positions in the (x, y) plane, where 
is maximum for a given z-value. 

From a computational point of view one have to set /3q 
lower than the solar values (IC* - 10^) to reduce the com- 
putational time-scale to a reasonable value. As a compro- 
mise I set the initial field strength so that /3o = 100 yield- 
ing Cm/gk = 100 with the present convection model. 
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